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Abstract 

We review several most recent prompt-photon measurements at HERA and the 
' (— I 'I Tevatron and discuss their implication for future measurements at the LHC. A 

comparison to Monte Carlo models, as well as to NLO QCD predictions based 
^ • on the standard DGLAP and the A;^ -factorization approaches is discussed. Ef- 

O I fects from renormalization and factorization scale uncertainties, as well as un- 

certainties on the gluon density distribution inside a proton are discussed. 



>- ■ 1 Introduction 

^ I Events with an isolated photon are an important tool to study hard interaction processes since such 

OO ■ photons emerge without the hadronization phase. In particular, final states of ep and pp collisions with 

a prompt photon together with a jet are more directly sensitive to the underlying QCD process than 
^ I inclusive prompt photon measurements. 

00 i The results on prompt-photon production provided by 

O ' HERA are important for the interpretation of the LHC data. Un- g 

like pp collisions, ep collisions involve a point-like incoming lep- 
^ I ton, which leads to some simplification in the description of the 

' prompt-photon production in ep compared to pp. At HERA, the 

quark content of the proton is probed through the elastic scatter- 
ing of a photon by a quark, — > -yq (see Fig. [T]). Such QED 
events are significantly simpler than lowest-order compton-like 
Q9 ~^ 11 events which are dominant in pp collisions (see Fig. [2l 
left figure). The latter process has direct sensitivity to the strong 
coupling constant and requires much better understanding of the 

gluon structure function inside both incoming protons than for the Fig. 1: Lowest-order diagram (Compton 
lowest-order diagram in ep collisions. scattering) for 7+jet events in ep collisions 

Despite the difference between ep and pp collisions con- ■ 
cerning certain lowest-order diagrams, a large class of partonic 

contributions are similar between ep and pp collisions, due to the hadronic nature of the resolved photon. 
In particular, a contribution to prompt-photon events from the gq — > qy process in photoproduction, 
in which the photon displays its hadronic structure [1-4], leads to significant sensitivity to the gluon 
structure function as is the case in pp collisions (see Fig. |2l the two figures on the right). Therefore, 




analysis of HERA data can make a bridge between a better understood ep case and the less understood 
pp case, since apart from the convolution with different structure functions, photoproduction diagrams 
ep collisions involving a resolved photon are essentially the same as diagrams in pp collisions. 



Fig. 2: The dominant diagram for prompt-photon events in pp collisions (left figure) compared to two resolved photon diagrams 
in ep photoproduction (see Section|2]for more details). 

Prompt-photon events in ep collisions can constrain both quark and gluon parton densities (PDFs). 
In addition, differences between coUinear factorization and kx factorization in the description of the un- 
derlying hard subprocesses can be studied in detail. This is important not only for a better understanding 
of QCD dynamics, but also has direct implications for searches of exotic physics at the LHC, in which 
prompt-photon production is the main background. A number of QCD predictions [2-6] can be con- 
fronted with the data and some of them will be described in more detail below. 

2 Photoproduction of prompt photons at NLO 

In the photoproduction ep scattering processes, the electron is scattered at small angles, emitting a quasi- 
real photon which scatters with the proton. The spectrum of these photons can be described by the 
Weizsacker- Williams approximation [7]. The photons will take part in the hard interaction either di- 
rectly, or through their "partonic" content, in which case a parton stemming from the resolved photon 
participates in the hard subprocess. Similarly, a photon in the final state can either originate directly from 
the hard interaction or from the fragmentation of a parton. Therefore, one can distinguish four categories 
of subprocesses: 1) direct direct, 2) direct fragmentation, 3) resolved direct and 4) resolved fragmen- 
tation. Examples of leading order diagrams of each class are shown in Fig. [3l Beyond leading order, 
this distinction becomes ambiguous. For example, the NLO corrections to the direct part involve final 
state collinear quark-photon pairs which lead to divergences which are absorbed into the fragmentation 
function, such that only the sum of these contributions has a physical meaning. The complete NLO cor- 
rections to all four parts have been calculated in [4] for inclusive prompt photons and in [8] for photon 
plus jet final states. A public program EPHOX, written as a partonic event generator, is available from [9]. 
The NLO corrections to the direct-direct part also have been calculated in [3, 10] for the inclusive and 
photon plus jet final state. 

The 7-p scattering processes are of special interest since they are sensitive to both the partonic 
structure of the photon as well as of the proton. They offer the possibility to constrain the (presently 
poorly known) gluon distributions in the photon, since in a certain kinematic region the subprocess 
qg — > 75, where the gluon is stemming from a resolved photon, dominates [5]. 

Working within the framework of collinear factorization, i.e. assuming that the transverse mo- 
menta of the partons within the proton can be neglected and other non-perturbative effects can be fac- 
torized from the hard interaction at high momentum transfers, the cross section for ep —>■ can 
symbolically be written as a convolution of the parton densities for the incident particles (respectively 
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fragmentation function for an outgoing parton fragmenting into a photon) with the partonic cross section 
a: 



da'P-*^^{Pp,Pe,Pj) = Y^ j dxej dxpj dzF,/,{xe,M)Fi,/p{xp,Mp)D^/,{z,MF) 

a,b,c 

da^^^'^ {XpPp, XePe, P^/z, fi, M, Mp, Mf) , (1) 

where M, Mp are the initial state factorization scales, Mp the final state factorization scale, /i the renor- 
malization scale and a, b, c run over parton types. In the NLO calculations shown in Fig. HI all these 
scales are set equal to pj, and varied simultaneously. The functions Fiy/p{xp, Mp) are the parton distribu- 
tion functions in the proton, obeying DGLAP evolution. Note that including initial state radiation at NLO 
in the partonic calculation means that the partons taking part in the hard interaction can pick up a nonzero 
transverse momentum. In certain cases, this additional "fcy-kick" seems to be sufficient to describe the 
data well. For example, a study of the effective transverse momentum {kx) of partons in the proton 
has been made by ZEUS [11]. Comparing the shapes of normalized distributions for (A;^) -sensitive 
observables to an NLO calculation, it was found that the data agree well with NLO QCD without extra 
intrinsic (kx) [8]. 





rcsolvcd-dircct resolved-fragmentation 
Fig. 3: Examples of contributing subprocesses at leading order to each of the four categories in ep collisions. 



The "resolved" contributions are characterized by a resolved photon in the initial state where a 
parton stemming from the photon instead of the photon itself participates in the hard subprocess. In 
these cases, Fa/e{xe, M) is given by a convolution of the Weizsacker-WiUiams spectrum fj{y) with the 
parton distributions in the photon: 

Fa/e{xe,M)= f dy dx^ f^{y) F^j^ix^, M) 5{x^y - x^) . (2) 



The cases with "direct" attributed to the initial state photon correspond to a = 7, so F^/^ = (5(1— x^) and 
F„/e in eq. ^ collapses to the Weizsacker- Williams spectrum. The cases "direct-direct" and "resolved- 
direct" correspond to c = 7, so D^i^{z, Mp) = ^^^^(l — z) in ([T]), i.e. the prompt photon is produced 
directly in the hard subprocess and not from the fragmentation of a hard parton. 

If additional jets are measured, eq. ([T]) also contains a jet function, which defines the clustering 
of the final state partons other than the photon into jets. Prompt photon production in association with 
a jet offers more possibilities to probe the underlying parton dynamics. It allows for the definition of 
observables that provide information about the longitudinal momentum fractions x'^,x^ carried by the 
particles taking part in the hard interaction. The partonic x^, x^ are not observable, but one can define 
the observables 
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which, for direct photons in the final state, coincide with the partonic x"* , x^ at leading order. Unique 
to photoproduction processes is the possibility to "switch on/off the resolved photon by suppress- 
ing/enhancing large x'^ . A?, x'^ = 1 corresponds to direct photons in the initial state, one can obtain 
resolved photon enriched data samples by placing a cut x^^^ < 0.9. Another possibility to enhance or 
suppress the resolved photon component is to place cuts on pT and rapidity. From eq. (O one can easily 
see that x^^^ is small at low values and large rapidities. Small x'^'-enriched data samples could be 
used to further constrain the parton distributions in the real photon, in particular the gluon distribution, 
as investigated e.g. in [5]. Similarly, one can suppress the contribution from the resolved photon to probe 
the proton at small x^ by direct 7-p interactions [5]. 

In order to single out the prompt photon events from the background of secondary photons pro- 
duced by the decays of Ught mesons, isolation cuts have to be imposed on the photon signals in the 
experiment. A widely used isolation criterion is the following: A photon is isolated if, inside a cone 
centered around the photon direction in the rapidity and azimuthal angle plane, the amount of hadronic 
transverse energy Elj?''^ deposited is smaller than some value i?T,max '■ 

for (7?-r/^)2 + {(t>-4>^f<R, 

TTihad ^ TP ^ ' 

— -CT.max • 

HERA experiments mostly used E'r.max = ^Pr with e = 0.1 and i? = 1. Isolation not only reduces the 
background from secondary photons, but also substantially reduces the contribution from the fragmen- 
tation of hard partons into high-pT photons. When comparing the result of partonic calculation to data, 
photon isolation is a delicate issue. For example, a part of the hadronic energy measured in the cone 
may come from the underlying event; therefore even the direct contribution can be cut by the isolation 
condition if the latter is too stringent. 



3 fcr -factorization approach 

A complementary description is offered by the /cr-factorization approach [12], which relies on parton 
distribution functions where the /cr-dependence has not been integrated out. 

In the framework of /c^— factorization approach the treatment of fcy -enhancement in the inclusive 
prompt photon suggests a possible modification of the above simple kx smearing picture. In this ap- 
proach the transverse momentum of incoming partons is generated in the course of non-coUinear parton 
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Fig. 4: The differential 7+jet cross sections as functions of Et and 77 of the prompt photon and the jet. The data are com- 
pared to QCD calculations and Monte Carlo models as described in the text. The shaded bands correspond to a typical scale 
uncertainty which was obtained by changing the renormalization and factorization scales simultaneously by a factor of 0.5 and 
2 respectively. 

evolution under control of relevant evolution equations. In the papers [6, 13] the Kimber-Martin-Ryskin 
(KMR) formalism [14] was applied to study the role of the perturbative components of partonic krp in 
describing of the observed Et spectrum at HERA and Tevatron. The proper off-shell expressions for 
the matrix elements of the partonic subprocesses and the KMR-constructed unintegrated parton densities 
obtained independently were used in [13] to analyze the Tevatron data. 

4 Comparison with HERA results 

Recently published [15] ZEUS differential cross sections as functions of Et and -q for the prompt- 
photon candidates and for the accompanying jets have revealed some difference with both Monte Carlo 
predictions and the next-to-leading order (NLO) calculations based on the collinear factorization and the 
DGLAP formalism [3,4], as shown in Fig. |4] The data are compared to QCD calculations performed 
by Krawczyk and Zembrzuski (KZ) [3], by Fontannaz, Quillet and Heinrich (FGH) [4], by A. Lipatov 
and N. Zotov (LZ) [6] and and Pythia 6.4 [16] and Herwig 6.5 [17] Monte Carlo models. The MC 
differential cross sections do not rise as steeply at low Ej, as do the data. It should be pointed out 
that no intrinsic transverse momentum of the initial-state partons in the proton was assumed for these 
calculations. The QCD calculation [6] based on the -factorization [12] and the Kimber-Martin-Ryskin 
(KMR) prescription [14] for unintegrated quark and gluon densities, gives the best description of the Et 
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Fig. 5: The differential cross section da/dEr and da /drf as functions of iSj and rf' of the inclusive prompt photon photo- 
production calculated at —0.7 < if < 0.9 and 0.2 < y < 0.9. The data are compared to two different NLO calculations. 
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Fig. 6: Same as in Fig. [5] but for 7+jet events with the additional jet cuts: — 1 < 77^"^* < 2.3 and E'^^ > 4.5 GeV. 



In the photon-rapidity distribution of Fig. IH the data lying above the NLO theory prediction at 
low values of rj"^ could be explained by the fact that in this region, x^j^^ is small, as can be seen from 
Eq. Q, and therefore A^^-effects may be important. On the other hand, this is not corroborated by the 
jet rapidity distribution, which has a problem at high rf°^, corresponding to small x^f^^. Indeed, a direct 
measurement [15] of x^^^ shows that the differences with NLO are mainly at low values of the x^i^g 
distribution. In this region, resolved photon events dominate, which may indicate that resolved photon 
remnants could have lead to an increase in the number of jets which have passed the experimental cuts, 
while these events are not accounted for in the partonic calculation. 

The inclusive prompt photon data [18, 19] lie above the NLO theory prediction in the whole ra- 
pidity range, except for the bin of largest rapidity, where the agreement is good after hadronization 



corrections, see Fig. [51 

Interestingly, ZEUS investigated what happens if the minimum transverse energy of the prompt 
photon is increased to 7 GeV, and found that in this case, the NLO calculations are in good agreement 
[15] , which suggests that non-perturbative effects may produce the discrepancy. See [20] and references 
therein for more details. 

The HI experimental data in photoproduction [18] are shown in Figs. [5]and[6l Both inclusive 
and 7+jet cross sections were compared to the FGH NLO calculations after hadronization corrections. 
The HI data [18] referred to the kinematic region defined by 5 < E'J < 10 GeV, -1 < -rf < 0.9 
and 0.2 < y < 0.7, which is rather similar to the ZEUS measurement shown in Fig. [4] Similar to the 
ZEUS case, MC predictions were found to underestimate the HI cross sections, while NLO QCD gives 
a much better description. After taking into account hadronization and multiple interaction effects, NLO 
calculations predict somewhat smaller cross sections compared the measurements [18]. 





Fig. 7: The differential cross section da/dEr and da /drf as functions of iSj and rf of tfie inclusive prompt photon photo- 
production calculated at —0.7 < rj'' < 0.9 and 0.2 < y < 0.9. The data are compared to the fcr— factorization calculations. 
The bands correspond to a typical renormalization scale uncertainty which was obtained by changing /ir by a factor of 0.5 and 
2. 





-0.8 -0.4 0.0 0.4 



Fig. 8: Same as Fig.|71 but for 7-t-jet events with the additional jet cuts: -1 < ry^''* < 2.3 and E'°^ > 4.5 GeV. 



The HI experimental data in photoproduction [18] were also compared to the /cj-factorization 
approach [6]. Comparison with the factorization approach indicates somewhat better agreement, as 
shown in Fig. [7] (see [6] for details). One can see that the measured distributions are reasonably well 
described except the moderate Ej, region and in the pseudorapidity region —0.4 < if < 0.9 only. 
For —1 < -rf < —0.4 the A;^— factorization predictions are mostly below the experimental points. The 
discrepancy between data and theory at negative r/"^ is found to be relatively strong at low values of the 
initial photon fractional momentum y. The effect of scale variations in transverse energy distributions is 
rather large: the relative difference between results for fi = Ej, and results for fi = Ej,/2 or fj, = 2Ej, is 
about 15% within the fc^-factorization approach, which is due to missing higher order corrections. The 
scale dependence of the NLO QCD calculations in the coUinear factorization approach is below the 10% 
level. 

The individual contributions from the direct and resolved production mechanisms to the total cross 
section in the fc-p -factorization approach is about 47% and 53%, respectively. The contributions of single 
resolved processes 

q{ki) + g{k2) -f{pj) + q{p'), 

g{ki) + q{k2) -f{p^) + q{p'), 

q{ki) + q{k2) ^ l{p^) + g{p'). 

account for 80%, 14% and 6% respectively. 

The transverse energy Ej, and pseudorapidity if distributions for 7+jet events measured by HI 
are compared to the fcy -factorization predictions in Fig.[8](see also Ref. [6]). In contrast to the inclusive 
case, one can see that the /c^-factorization predictions are consistent with the data in most bins, although 
some discrepancies are present. The theoretical results are lower than the experimental data at negative 
7]'^ and higher at positive r/'^. The scale variation as it was described above changes the estimated cross 
sections by about 10%. Note that such disagreement between predicted and measured cross sections has 
also been observed for the NLO QCD calculations in the coUinear factorization approach, see Fig.|6] 




Fig. 9: The differential cross section da / dx^ and da /dxp of prompt photon + jet production calculated at 5 < iJj < 10 GeV 
and 0.2 < y < 0.7 with an additional jet requirement -1 < 7)'"^ < 2.3 and Ei^^ > 4.5 GeV. 



Figure |9] shows the x^f^^ and x^f^^ distributions (see Eq. O measured by HI. One can see that kx 
factorization predictions reasonably well agree with the experimental data. The NLO calculations [3,4] 
without corrections for hadronization and multiple interactions give similar results. 



The HI Collaboration [21] also has performed 7+jet measurements in DIS for > A GeV^. The 
NLO calculations [22], which are only available for 7+jet final state, failed to describe normalization 
of the cross sections, although the agreement in shape was found to be reasonable (Fig. [TOl ). No kx 
factorization prediction available for DIS. 
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Fig. 10: The differential 7+jet cross sections as functions of Et and 77 of the prompt photon in DIS. The data are compared to 
LO and NLO calculations. 



In summary, some differences with NLO QCD were observed in both photoproduction and DIS. 
Differences at low Pj, can be due to the treatment of the fragmentation contribution in NLO calculations. 
Further, it would be interesting to see the effect of calculations beyond NLO QCD. The approach based 
on the kx factorization has better agreement with the data, but such calculations have larger theoretical 
uncertainties. 



5 Comparison with Tevatron results 

Isolated photons in pp collisions at Tevatron have been measured recently by the CDF [23, 24] and 
DO [25-28] Collaborations. 

Measurements of pp — > 7 + jet + X for 30 GeV < Pj. — 300 GeV have very recently been 
published by DO [28]. The comparison to theory is done separately for different regions in rapidity of 
the photon and the jet. The NLO partonic Monte Carlo program JETPHOX [9, 29] was used to compare 
the data to theory at next-to-leading order. It was shown that the NLO calculations are not sufficient 
to describe the shape of Pj, distributions in different rapidity regions, as can be seen in Figure [TT] At 
present, the comparison with the /cy -factorization prediction is in progress. 

Differences with the collinear factorization approach have been seen previously as well. Both 
CDF [23] and DO [26] cross sections were found to be aboveQ NLO predictions at low P^. However, 
RHIC has also measured prompt photon production in pp collisions at ^/s = 200 GeV and found good 
agreement with NLO theory in the collinear factorization approach [20, 30]. 

The same data were compared to the kx factorization approach in [13]. Figures[T2]and [T3]show the 
CDF [23] and DO [26] measurements for the da / dEj,drj'^' cross sections calculated at ^/s = 630 and 1800 
GeV in central and forward kinematic regions together with the kx factorization predictions. One can 
see that theoretical predictions agree with the experimental data within the scale uncertainties. However, 

'For DO, the difference was mainly concentrated in the central rapidity region. 
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Fig. 11: The ratios of the triple-differential cross section measured by DO compared to to the NLO QCD prediction using 
JETPHOX. See details in Ref. [28]. 




Fig. 12: The double differential cross section da / dE'^drp of inclusive prompt photon production at -y/s = 630 GeV and 
1 77''! < 0.9 (left plot) and 1.6 < \rp\ < 2.5 (right panel). The solid line corresponds to the default scale ^ = i5j of the kr 
factorization predictions, whereas upper and lower dashed lines correspond to the fj, = _Ej/2 and /i = 2_EJ. 



the results of the calculation with the default scale tend to underestimate the data in the central kinematic 
region and agree with the DO data in the forward 'q'^ region. The coUinear NLO QCD calculations give 
a similar description of the data: generally there is a residual negative slope in the ratio of the data over 
the prediction as a function of Ej,. The scale dependence of the factorization results is rather large 
(20 — 30%), due to the fact that these are leading order calculations. 
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Fig. 13: The double differential cross section da / dE'^dif of inclusive prompt photon production at ^/s — 1800 GeV and 
1 77^ I < . 9 (left plot) and 1.6 < l?;^] < 2.5 (right plot) . The solid line corresponds to the default scale ^ — , whereas upper 
and lower dashed lines correspond to the /i = iJj /2 and ^ — 2E^ for the kr factorization calculations. 



The double differential cross sections d(7/dEj.dvf are usually the most difficult observables to 
describe using QCD predictions. Yet, as it can be seen from Fig. [T3j the /c-r— factorization predictions 
agree well with DO [25] and CDF [23] data both in shape and normalization. There are only rather 
small overestimations of the data at low E'J values in Figs. [13] in the forward region. Again, the scale 
dependence of our calculations is about 20-30%. The theoretical uncertainties of the collinear NLO 
predictions are smaller (about 6% [25]), which is to be expected as inclusion of higher order terms 
reduces the scale uncertainty. 

One can conclude that the results of calculations in the ^7-— factorization approach in general agree 
well with Tevatron experimental data, within a large scale uncertainty. 



-NLOCTEQ6.I 14 TeV 



NL0CTEQ6.1 1.96 TeV 



6 Prompt photons at LHC 

The direct photon production at LHC has significantly higher 
cross sections compared to the ones measured at Tevatron and 
HERA. The prompt-photon cross section at LHC is more than a 
factor of hundred higher than that at Tevatron and a factor of 10^ 
larger than that for photoproduction at HERA, assuming a similar 
kinematic range (| if |< 2), as shown in Fig. [141 This will al- 
low to explore the TeV energy scale already in a few years of data 
taking. 

Figure [15] shows the comparison between PYTHIA and 
HERWIG Monte Carlo models and JETPHOX LO and NLO cal- 
culations. The cross sections for 7 + jet events were calculated 
for I |< 2, > 100 GeV and > 105 GeV. The cuts 
on the transverse momenta are asymmetric to avoid instabilities 
in the NLO calculations. An isolation requirement i?^ > 0.9 E^^ 
was imposed, where Ej^^ is the total energy of the jet which con- 
tains prompt photon. Jets were reconstructed with the longitudinally-invariant kx algorithm in inclusive 
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Fig. 14: The PJ. cross section for 7+jet 
events predicted by NLO QCD for the Teva- 
tron and the LHC kinematic range. 



mode [31]. 

The NLO QCD calculation is 30^0% higher than that predicted by PYTHIA. On the other hand, 
PYTHIA is 20% above HERWIG. It is interesting to observe that the level of discrepancy between 
PYTHIA and HERWIG is about the same as that observed at HERA at much lower transverse momenta 
(for example see Fig. |4ll. However, there is no significant difference between NLO and PYTHIA at 
> 10 GeV for ep, while at the LHC energy range the difference between NLO and PYTHIA is rather 
significant. Certainly, the overall normalization of Monte Carlo programs like PYTHIA or HERWIG has 
to be adjusted, as these programs cannot account for contributions from loop corrections at higher orders. 




Fig. 15: Comparisons of Monte Carlo models with LO and NLO calculations as implemented in JETPHOX. 



Scale uncertainties were estimated by changing the renormalization and factorization scales in the 
range 0.5 < < 2. The relative difference between predicted cross sections is shown in Fig. [T6l 

To make quantitative statements on scale uncertainties with the present level of statistical errors in calcu- 
lations using JETPHOX, a linear fit was performed to determine the trend of the relative differences with 
increase of Pj^. As it can be seen, the scale uncertainty is about 10% and slowly increases with Pji . 

To estimate the uncertainty associated with the gluon density, the calculations have been performed 
using two CTEQ6.1M sets (15 and 30) which correspond to two extremes in the gluon density at large 
X [32]. Fig. [T7] shows the relative difference between those two sets as a function of i-^. It is seen 
that the gluon uncertainty is almost a factor of two larger compared to the scale uncertainty estimated 
above. No statistically significant difference has been observed between the cross sections calculated 
using CTEQ6. IM and MRST04. This is not totally surprising as both sets have similar input data for the 
global fit analysis. 

The predictions for the kx factorization approach were obtained for a wider pseudorapidity range, 
for both central and forward pseudo-rapidities r/'^. As a representative example, we will define the central 
and forward kinematic regions by the requirements \7]'^\ < 2.5 and 2.5 < \rj'^\ < 4, respectively. The 
transverse energy E'J distributions of the inclusive prompt photon production in different r/"^ ranges at 
= 14 TeV are shown in Figs. [18] One can see that variation in scale /x changes the estimated cross 




Fig. 16: Relative difference between the cross section estimated 
with fi = 0.5 (P^") and ^ = 2 (P|°™") as a function of Pt 
for gamma and jet. The line represents a linear fit. 



Fig. 17: Relative difference between the cross section estimated 
with CTEQ6.1M set=-15 (P^^) and set =15 (P|°'"'") as a func- 
tion of Pt for gamma and jet. The line represents a linear fit. 



sections by about 20-30%. However, as it was already discussed above, there are additional theoretical 
uncertainties due to the non-coUinear paiton evolution, and these uncertainties are not well studied up to 
this time. Also the extrapolation of the available parton distribution to the region of lower x is a special 
problem at the LHC energies. In particular, one of the problem is connected with the correct treatment 
of saturation effects in small x region. Therefore, more work needs to be done until these uncertainties 
will be reduced. 

Thus, the calculation based on the factorization approach shows a larger scale uncertainty 
compared to the coUinear factorization approach: for Pj, ~ 100 GeV, the overall uncertainty for the NLO 
calculations is expected to be around 10%, while it reaches 20-30% for the /ct -factorization calculations 
for the same Pj, range, due to the fact that the latter are at leading order in as- As the residual scale 
dependence of missing higher order terms resides in logarithms involving ratios of Pj' and scales jj?, the 
effect becomes more dramatic at the LHC energy. 

7 Summary 

In this review, we have attempted to summarize recent progress in the description of prompt photon 
production at HERA, the Tevatron and the LHC. At HERA, some differences with NLO were observed 
in both photoproduction and DIS. The deficiencies at low Pj, values may indicate that non-perturbative 
effects at small Pj, play a non-negligible role. Also, one should expect that adding high-order correc- 
tions to the collinear-factorization approach should improve the description. Similar conclusions can be 
drawn for the Tevatron data which, as in the HERA case, has differences with NLO in the lowest Pj, 
region. Recently, significant differences with NLO were observed by the Tevatron for the shapes of Pj, 
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Fig. 18: Left plot: The kr factorization predictions for differential cross sections da / dEi^drp at = 14 TeV GeV and 
\rf' I < 2.5 (left plot); at 2.5 < rp < 4.0 (right plot). The solid line corresponds to the default scale /^ = E'^, whereas upper 
and lower dashed lines correspond to the yu = E"^ /2 and /j = 2i?J. 

distributions differential in rf. On the other hand, RHIC observes good agreement with NLO QCD. 
Considering the fact that RHIC uses a photon isolation method which is different from the usual cone 
isolation, the differences mentioned above may also have to do with isolation criteria acting differently 
in a partonic calculation than in the full hadronic environment of the experiment. 

An alternative approach based on the Ut factorization generally improves the description of the 
HERA and the Tevatron data, but it has larger theoretical uncertainties. As for NLO, high-order correc- 
tions to the /cT-factorization approach should improve the description of the data. The apphcabiUty of 
the kx factorization to the LHC data will be tested with the arrival of the first LHC data, but it is already 
evident that significant theoretical uncertainties are expected for the description of prompt-photon cross 
sections at LHC. Using the the coUinear factorization approach, uncertainties of NLO calculations are 
expected to be 10-20% at about 1 TeV photon transverse momenta, and significanfly larger for the kr- 
factorization calculations. These uncertainties have to be reduced in the future for detailed comparison 
of the LHC data with the QCD predictions. 

In all cases, Monte Carlo predictions fail to describe prompt-photon cross sections, both in shape 
and normaUzation. Generally, HERWIG is significantly below PYTHIA. This could have a direct impact 
on the future LHC measurements, in particular for exotic searches which often rely on Monte Carlo 
predictions for estimations of rates for background events. 
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